function meanerr(time,y,R,sig,y0)

load('orthbase');

exactmean=zeros(1,100);
time(1)=10^-15;

for ii=1:100
    exactmean(ii)=y0(1)+quadl(@(s) R(s)-1/2*sig(s).*sig(s),0,time(ii));
end
hold on
plot(time,exactmean);
plot(time,y(:,1),'r+');
title('Mean')
legend('exact mean','wce mean')
